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Abstract 

The doping of graphene to tune its electronic structure is essential for its further use in carbon 
based electronics. Adapting strategies from classical silicon based semiconductor technology, we use 
the incorporation of heteroatoms in the 2D graphene network as a straightforward way to achieve 
this goal. Here, we report on the synthesis of boron-doped graphene on Ni(lll) in a chemical vapor 
deposition process of triethylborane on the one hand and by segregation of boron from the bulk 
on the other hand. The chemical environment of boron was determined by x-ray photoelectron 
spectroscopy and angle resolved photoelectron spectroscopy was used to analyze the impact on 
the band structure. Doping with boron leads to a shift of the graphene bands to lower binding 
energies. The shift depends on the doping concentration and for a doping level of 0.3 ML a shift 
of up to 1.2 eV is observed. The experimental results are in agreement with density-functional 
calculations. Furthermore, our calculations suggest that doping with boron leads to graphene 
preferentially adsorbed in the top-fcc geometry, since the boron atoms in the graphene lattice are 
then adsorbed at substrate fcc-hollow sites. The smaller adsorption distance of boron compared to 
carbon leads to a bending of the graphene sheet in the vicinity of the boron atoms. By comparing 
calculations of doped and undoped graphene on Ni(lll), as well as the respective free-standing 
cases, we are able to distinguish between the effects that doping and adsorption have on the band 
structure of graphene. Both, doping and bonding to the surface, result in opposing shifts on the 
graphene bands. 
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ii I. INTRODUCTION 



12 The concept of doping was initially introduced in classical semiconductor technology, 

13 but was also adapted for carbon allotropes^l Doping of carbon material with nitrogen 
u and boron atoms generated intense interest due to the possibility of tailoring the physical 

15 properties, i.e., electronic and transport properties'^ but also chemical properties like the 

16 ability to adsorb lithium for capacitors^ or hydrogen storage capabilities^. In particular, 

17 these concepts were discussed recently for the two dimensional carbon allotrope graphene. 
is The doping of graphene with nitrogen was achieved from nitrogen doped precursors or from 

19 post growth treatment procedures'^. The resulting new material showed promising first 

20 results towards applications in the field of electrochemical sensin gECSEI lithium batteries 131 , 

21 material in p-n junctions^, and fuel cells^. Nitrogen doping of graphene showed a fun- 

22 damental dependence on the geometry/site of the dopant. Besides the expected n type 

23 behavior for substitutional doped graphene, p type doping was observed for a dopant ge- 

24 ometry, where nitrogen is introduced next to a carbon vacancy, forming pyridine-like units 

25 within the graphene lattice^EH. The changes in the band structure upon nitrogen-doping 

26 showed, nevertheless, the possibility of tuning the graphene band structure. The thermal 

27 stability of the resulting nitrogen-modified graphene sheets is similar to that of graphene. 

28 Based on these results the incorporation of boron was also considered and first results from 

29 boron-doped graphene were already presente d 22 1 23 1 . In both casese exfoliated graphene was 

30 used, implying the already known challenges for the quality of the graphene layers and for 

31 large scale production. Also first results from the incorporation boron in grapheneoxide are 

32 reported.^! 

33 Herein, we report the production of single layer boron-doped graphene on a Ni(lll) sur- 

34 face by chemical vapor deposition (CVD) of the boron-containing precursor triethylborane 

35 (TEB) and from segregation. The two processes used yield tunable concentrations of boron 

36 in a highly ordered graphene layer. The boron-doped graphene is characterized with angle 

37 resolved photoemission spectroscopy (ARPES) and x-ray photoelectron spectroscopy (XPS). 

38 The experimental results are analyzed and discussed by comparison with density-functional 

39 theory (DFT) calculations of undoped an boron-doped graphene on either Ni(lll) or free- 

40 standing in vacuum, respectively. Our data show a strong shift of the graphene bands of 

41 up to 1.2 eV with respect to the Fermi level, depending on the boron concentration. Dur- 
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42 ing the doping process the overall band structure is retained. The comparison to the DFT 

43 calculations shows excellent agreement with the experiments. Our calculations also suggest 

44 that doping with boron leads to graphene being adsorbed in the top-fcc geometry, because 

45 the boron atoms in the graphene lattice preferentially adsorb in the substrate fcc-hollow 

46 sites. The smaller adsorption distance of boron atoms within graphene compared to carbon 

47 atoms leads to a bending of the graphene sheet in the vicinity of the boron atoms. By com- 

48 paring free-standing and/or undoped graphene with the adsorbed boron-doped graphene 

49 we are able to discuss the effects that doping and adsorption have on the band structure 

50 separately. Both modifications result in opposing shifts on the graphene bands. For high 

51 dopant concentrations this results in a Dirac point above the Fermi level. 

52 II. EXPERIMENTAL SETUP AND COMPUTATIONAL DETAILS 

53 The experiments were conducted at the third generation synchrotron source BESSY II 

54 in Berlin, Germany. Boron-doped graphene on Ni(lll) was produced by CVD using TEB 

55 as precursor at temperatures between 600 and 950 K, leading to boron concentrations of 

56 typically 0.15 up to 0.35 ML after an exposure of 1800 L. Concentrations below 0.15 ML 

57 were prepared by segregating boron from the bulk, while exposing the nickel crystal to 

58 propene at 10" 6 mbar at 900 K until saturation of the carbon signal. The boron is dissolved 

59 in the bulk by exposure to TEB and subsequent annealing to temperatures of 1100 K. The 

60 ARPES measurements were carried out at beamline U 56/2 PGM 2 using a Phoibos 100 
ei analyzer^, while a transportable set up was used for the XPS measurements at beamline 

62 U 49/2 PGM 1^. All XP spectra are taken with a photon energy of 380 eV and an overall 

63 resolution of 200 meV at normal emission. The boron concentration was calibrated from the 

64 adsorption of TEB at 130 K: we compared the TEB carbon intensity to the known intensity 

65 of a saturated benzene layer on Ni(lll) at 200 KPh From this we calibrated the boron 
ee coverage from the known boron content in TEB. The ARPES data was recorded using a 

67 photon energy of 70 eV. The DFT calculations were carried out using the VASP program 

68 package^, employing a plane wave basis and the Projector- Augmented- Wave (PAW) method 

69 for treating the core electrons^. The Perdew-Burke-Ernzerhof (PBE) 29 functional, including 

70 a correction to take dispersive forces into account!^, was applied. This theoretical setup was 

71 shown to accurately describe the interaction of graphene with metal surfaces^®!. The 
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72 graphene /metal systems were modeled by six layer slabs of nickel. During optimization, 

73 the topmost three nickel layers were relaxed to model the surface, while the bottom three 

74 layers were fixed at the calculated nickel bulk positions. The unit cell of the adsorbed 

75 graphene containing two carbon atoms was, in accordance with experiments, chosen to be 

76 commensurable with the nickel surface unit cell containing one nickel atom. Different boron 

77 coverages were simulated by replacing one carbon atom by a boron atom in the graphene 

78 sheet of (2x2), (3x3), and (4x4) graphene unit cells, leading to dopant concentrations 

79 with respect to the nickel atoms in the first substrate layer of 0.25, 0.11, and 0.063 ML, 
so respectively. In the (2x2) cells 13 x 13 x 1 Monkhorst-Pack 34 k-point grids and respective 
si smaller ones for the larger unit cells were used. They ensured a convergence with respect 

82 to total energies and band energies within ±0.01 eV. For further details on the theoretical 

83 setup se< 



84 III. RESULTS AND DISCUSSION 



85 A. Preparation, Doping Geometry, and Doping Level 
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FIG. 1. High resolution XP spectra of (a) the C Is and (b) the B Is region of boron-doped graphene 
with different doping levels. The inset in (a) shows the shift of the graphene C Is peak with the 
doping concentration. 



ae Figs. [T] (a) and (b) show typical C Is and B Is XP spectra of boron-doped graphene 

87 layers grown on Ni(lll). In the C Is region two main peaks are observed, the graphene 

88 peak at ~285.0 eV and a second one shifted to lower binding energies by about ~1.6 eV. 

89 For further analysis, we calculated core level shifts (CLS) of the C Is and the B Is states 



90 according to the final state approximation^, for different possible arrangements (Tab. [T]). 

91 From comparison of the calculated and the experimentally observed shifts in the C Is spectra 

92 the peak at lower binding energies is attributed to the formation of C 2 CB, i.e., substitutional 

93 boron (theoretical shift of 1.42 eV). The shift towards lower binding energies is caused by 

94 the increased charge density on the carbon atoms due to their larger electronegativity in 

95 comparison to the adjacent boron atom. This is also in line with a Bader analysis^ carried 

96 out for the calculated systems that yields a charge transfer of 1.8 e from each boron atom to 

97 the neighboring carbon atoms. The formation of a new peak in the C Is region is in line with 

98 the results obtained for nitrogen doped graphene, where a new signal was observed at about 

99 0.7 eV higher binding energies compared to the bridge-top graphene peak^. These shifts 

100 are expected considering the electronegativities of carbon, nitrogen, and boron. Besides the 

101 formation of this additional peak in the C Is spectrum, due to the bonds of carbon to the 

102 electropositive boron atoms, we additionally observe a shift of the main peak from 285 to 

103 284.4 eV with rising boron coverage (see inset of Fig. [T] (a)). The reason for this is that 

104 the boron dopants are influencing not only the C Is level of direct neighbors. A shift is 

105 also observed in the CLS analysis of carbon atoms in next neighbor spheres (Tab. [T]). This 
loe means that boron has also an influence on C3C carbon atoms resulting in a shift and a 

107 broadening of the C Is peak with increasing doping concentration. At even higher boron 

108 coverages a third peak at 282.2 eV is rising. The assignment of this smaller contribution is 

109 not unambiguously, but according to our CLS calculation we assign it to carbon bound to 
no two substituted boron atoms (shift of 3.16 eV in Tab. [I]) Therefore, we suggest that primarily 
in substitutional boron-doped graphene is formed, which is in line with the typical coordination 

112 sphere of boron, the applied growth mechanism, and is confirmed by our DFT calculations 

113 discussed in the following. 

114 In the B Is region there are two main contributions, located at 187.4 eV and a shoulder 

115 at 188.0 eV, respectively. They are attributed to boron bound to carbon, i.e., the species 
lie C2CB, and elemental boron, respectively. Please note that according to the calculated CLS 
117 the binding energy of the B Is electrons of boron bound to a different amount of carbon 
us atoms (C2CB, CCB2, and CB3) are very similar. This makes a clear assignment of the 

119 boron geometry from the B Is CLS difficult and, therefore, such conclusions were drawn 

120 from the data of the C Is region. From a quantitative analysis of the C Is and B Is data the 

121 combined carbon and boron coverage on the surface was calculated to be 2 ML, as expected 
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122 for a closed monolayer of graphene on Ni(lll). 



Species C Is (eV) B Is (eV) 

C 2 CB: < C-B 1.42 (0.95) -0.81 



CCB 2 : \ C-B 3.16 -0.06 



CB,: \ C-B 4.17 0.41 




C 2 C: \ C 2.86 



C 2 CB-B: < C-B 2.04 0.53 

TABLE I. Calculated C Is (relative to graphene in top-fcc adsorption) and B Is (relative to a 
single boron atom) core level shifts for different geometrical arrangements within the (2x2) cell of 
graphene adsorbed top-fcc on Ni(lll). Different arrangements in the graphene network of a carbon 
atom with one, two, and three boron neighbors (C2CB, CCB2, and CB3), a carbon atom next to 
a one atom vacancy (C2C), and a carbon atom adjacent to a boron-boron bond (C2CB-B) were 
calculated. The sketches display the relevant neighbors within the graphene lattice. The calculated 
C Is shift always refers to the central carbon atom in the shown sketch (the value in brakets in the 
C2CB case refers to the shift of the carbon atoms bond to the central carbon, i.e., carbon atoms 
C3C with three neighboring carbon atoms). A positive shift refers to a shift towards lower binding 
energies. 
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123 B. Influence of the Boron Doping on the Geometry 

124 Since graphene is known to adsorb on Ni(lll) in bridge-top and top-fcc adsorption geome- 

125 tries with almost identical binding energies^^, both arrangements were considered in the 

126 calculations of graphene and substitutional boron-doped graphene in all unit cells. During 

127 geometry optimization, an interesting shift of the adsorption geometry was observed. While 

128 for undoped graphene both adsorption geometries turned out to be local minima on the 

129 potential energy surface, the top-fcc geometry is found to be the most stable arrangement in 

130 all doped cases. The geometry optimizations of substitutional boron-doped graphene, which 

131 were started in bridge-top geometry, converged to the top-fcc geometry. The energy gain of 

132 top-fcc compared to bridge-top adsorption is estimated to be 0.3-0.4 eV per unit cell in all 

133 cases, by comparing the optimized structures to calculations without geometry optimization 

134 in the respective fixed bridge-top adsorption geometries. This is due to the fact that boron 

135 atoms prefer to be located over a fcc-hollow site, leading to the top-fcc arrangement for the 

136 whole graphene sheet, although the carbon atoms are known to slightly prefer the bridge-top 

137 over the top-fcc arrangement^. Note, however, that the adsorption geometry could not be 

138 determined in our experiments, since the graphene C Is and the broad C2CB C Is peak are 

139 superimposed and the latter dominates the splitting of the graphene C Is peak in the case of 

140 the top-fcc geometry. Nevertheless, this geometry is found most likely in our preparations, 

141 since the boron is not introduced after graphene growth, but is present during growth. In 

142 addition, the observed shift from bridge-top towards top-fcc adsorbed graphene layers might 

143 also be possible in other preparation techniques, starting from primarily bridge-top adsorbed 

144 graphene, since the energetic difference between bridge-top and top-fcc adsorbed graphene 

145 is known to be smalP^. The graphene band structure, however, is known to be rather in- 

146 sensitive to the two possible geometric arrangements, neglecting fine differences in the close 

147 vicinity of the Dirac point and thus the observed shifts are equal for both geometries. This 
we is confirmed by carrying out band structure calculations also of the unrelaxed bridge-top 

149 structures. The observed band shift of the optimized top-fcc structures is qualitatively as 

150 well as quantitatively reproduced. The differences of <0.06 eV are negligible and are due 

151 to the smaller distance of the graphene layer to the substrate in the case of the optimized 

152 top-fcc calculations. An additional effect observed during geometry optimization is that 

153 the boron atom is, in the doped cases, adsorbed closer to the surface than carbon atoms 
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154 in undoped graphene that are adsorbed in 2.12 A from the nickel substrate. In the (2x2) 

155 cell (0.25 ML) the adsorption distance of the boron atom is 2.02 A, i.e., the distance to 

156 the substrate is decreased by 0.1 A. The doping by boron atoms also affects the carbon 

157 atoms, which, in the (2x2) cell, exhibit a decreased adsorption distance of ~2.09 A. In both 
we larger unit cells the boron atom is adsorbed even closer to the nickel surface (1.99 A). This 
159 leads to a bending of the graphene layer: In the (4x4) cell the carbon atoms are adsorbed 
wo at distances between 2.08 and 2.14 A, i.e., the carbon atoms that are directly bound to the 

161 boron atom are adsorbed closer to the surface, while the carbon atoms that are located 

162 further away from the boron atom are adsorbed at slightly increased adsorption distances, 

163 compared to undoped graphene. 

164 C. Influence of Doping and Adsorption on the Graphene Band Structure 

165 It is known that the energetic position of graphene bands on metal surfaces depends 
we crucially on the adsorption distance®^, i.e., the strength of the adsorbate-substrate in- 

167 teractions. Therefore, the effects on the band structure due to the above discussed small 

168 change in the adsorption distance in the case of doping and due to the incorporation of 

169 boron into the graphene geometry, were studied. To that end we compared the boron- 

170 doped graphene layer in the (2x2) cell with and without taking into account the geometry 

171 changes due to boron-doping in the two adsorption distances of the optimized substitutional 

172 graphene (2.08 A) and the optimized undoped graphene (2.12 A). Both effects, the change 

173 in the adsorption distance and the geometry changes within the graphene sheet, influence 

174 the graphene band structure only negligible < 0.1 eV. This means, we can exclude the possi- 

175 bility that the observed shifting of carbon bands is originating from geometrical changes due 

176 to the boron incorporation. Instead, this shift has to originate from the different chemical 

177 and physical properties of the boron dopants compared to the carbon atoms in the graphene 

178 network, in particular the effect of the dopants on the electronic structure. 

179 In the following, we first discuss the effect of boron-doping as observed in our experiments, 
lso Secondly, based on the DFT band structure calculations, we analyze in details, on the one 

181 hand, the impact of doping free-standing or adsorbed graphene and, on the other hand, of 

182 the effect of adsorbing pristine or doped graphene on the Ni(lll) substrate. 

183 In Fig. [2] the ARPES data of boron-doped graphene for different dopant concentrations 
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FIG. 2. ARPES measurement of boron doped graphene with boron concentrations of (a) 0.28 ML, 
(b) 0.09 ML, and (c) 0.045 ML. 

iB4 are displayed. Fig. [2] (c) corresponds to a doping of 0.045 ML. The typical band dispersion of 

185 graphene on Ni(lll) is observed 12 * 32 * 41 * ^, but the energetic position of the ir band is shifted 

186 by 0.16 eV to smaller binding energies, in comparison to undoped graphene on nickel. The 

187 nickel bands are also visible, especially in the binding energy region between and 2 eV, 

188 where they couple to graphene states of appropriate symmetry^. The spectra in Figs. [2] (a) 

189 and (b), which show higher doping concentrations, reveal a shift of the bottom of the 7r band 

190 to subsequently smaller binding energies, depending on the doping level, up to a maximal 

191 shift of 1.2 eV in the case of a doping with 0.29 MLs of boron (Fig. [2] (a)). A quantitative 

192 analysis of this shift is given in Fig. [3j Interestingly, we find that the spectra at higher 

193 doping levels become more diffuse, due to the incorporation of boron. In the data this is 

194 observed as a lower contrast in the shifted band structure. This is attributed to the fact that 

195 boron atoms, besides from doping the layer, are also defects leading to new states localized 
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FIG. 3. Shift of the binding energy Eb with the boron content of the carbon it band at the T 
point of graphene on Ni(lll) as measured in experiment (black) and predicted by theory (red). 

196 around the original graphene bands. Due to the shift of the graphene bands to lower binding 

197 energies, we also observe hybridization of the 7r band at the M point with nickel d bands. 

198 Note that we did not find a change in the band dispersion or similar effects in our analysis 

199 of the ARPES data, i.e., the shift to lower binding energies is the only observed effect. 

200 The experimental band structures in our study of boron-doped graphene adsorbed on 

201 Ni(lll) are altered in two different ways compared to free-standing graphene. That is, on 

202 the one hand, the effects on graphene exerted by the nickel substrate and, on the other hand, 

203 the effects on graphene due to the dopant atoms. In the following, these effects are discussed 

204 separately by considering the calculated band structures for the case of maximal doping, 

205 i.e., the (2x2) unit cells, where the effects are most distinct. Four cases are considered, see 

206 Fig. [5j (a) free-standing, pristine graphene, (b) free-standing, boron-doped graphene, (c) 

207 pristine graphene adsorbed top-fcc on Ni(lll), (d) boron-doped graphene adsorbed top-fcc 

208 on Ni(lll). In Fig. [5] band structures along the path marked in green in the reciprocal cell 

209 corresponding to the calculated (2x2) unit cell are shown, see Fig. [4] (b) for the different 

210 involved reciprocal cells and paths. 

211 The paths marked in Fig. [4] (b) in blue, cyan, and magenta, in the reciprocal cell corre- 

212 sponding to a (lxl) graphene unit cell are backfolded on the path marked as green arrow 

213 of the smaller (2 x 2) reciprocal unit cell which underlies Fig. [5j Note, however, that the k 

214 points denoted in the legend of Fig. prefer to the original, large (lxl) reciprocal cell of 
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FIG. 4. The sketch in the middle displays the reciprocal cell to a (lxl) hexagonal graphene unit 
cell (dark gray) within its Brillouin zone (black) and the respective reciprocal cell to a (2x2) unit 
cell (light gray). The unit cell vectors and the positions of high symmetry points M, M', K, 
K', and T are drawn orange and green for the reciprocal (lxl) and (2x2) unit cell, respectively. 
The bands along the paths colored in blue, cyan, and magenta in the reciprocal (lxl) unit cell 
are backfolded on the path marked as green arrow in the reciprocal (2x2) unit cell. To the left 
the band structure of free-standing, pristine graphene in the reciprocal (2x2) unit cell along the 
green marked arrow is displayed with the color coding designating the corresponding paths in the 
reciprocal (lxl) unit cell the bands are backfolded from. To the right the band structure of free- 
standing pristine graphene in the reciprocal (lxl) unit cell along the orange arrow is displayed 
with the color coding relating the bands to the corresponding ones displayed in the reciprocal (2x2) 
unit cell. (The cyan bands are the mirror image to those displayed in the reciprocal (2x2) unit 
cell.) 



215 graphene. (The F, M, and M' points of the reciprocal cell corresponding to the (lxl) unit 

216 cell are backfolded on top of each other in the reciprocal cell corresponding to the (2x2) unit 

217 cell.) For pristine graphene the bands along the two magenta- marked paths are identical 

218 due to symmetry. In case of doping this symmetry is destroyed and the bands along these 

219 two paths become different. This shows up in additional bands in Fig. [5] in the doped cases 

220 (b) and (d) compared to the undoped cases (a) and (c). The experimental band structure of 

221 Fig. [2] is the one along the path marked as orange arrow in Fig. [4] (b) in the (1 x 1) reciprocal 

222 unit cell of graphene. 

223 Comparison of the band structures of pristine and doped free-standing graphene, Fig. [5] (a) 

224 and (b) shows that doping has two effects. Firstly, the shape of all bands is somewhat 
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FIG. 5. Band structures of free-standing pristine graphene (a), free-standing graphene doped 
substitutionally by boron (b), pristine graphene adsorbed top-fcc on Ni(lll) (c), and boron-doped 
graphene adsorbed top-fcc on Ni(lll) (d) calculated in (2x2) unit cells. Red, orange, and yellow 
circles represent carbon 2s, 2p x /2p y , and 2p z contributions to carbon s, a, and ir bands, with 
their radius being correlated to the manitude of the contribution to the bands (the radii of p z 
contributions are doubled for better visualization). Color coding is different for pristine graphene 
calculated in vacuum (a). Here, carbon 2s, 2p x /2p y , and 2p z contributions to bands that are 
backfolded from paths colored in Fig. [4] (b) in magenta in the reciprocal cell corresponding to 
the (lxl) unit cell are marked by cyan, purple, and blue circles, respectively. Nickel bands are 
represented by gray lines. 

225 changed leading, among other things, to an opening of the Dirac point and, as discussed 

226 above, to the appearance of additional bands due to the lowering of the symmetry upon 

227 doping, i.e., upon replacing one carbon atom in the (2x2) unit cell by boron. Secondly, 

228 doping leads to a general shift of all graphene valence bands towards lower binding energies 

229 by about 1.2-2.0 eV with respect to the Fermi energy. This shift can be interpreted as 

230 classical doping effect. In the limit of classical doping the band structure is assumed to not 

231 change at all but the number of electrons in the system is modified. If we neglect the change 

232 of the bands in our case, then the replacement of carbon by boron decreases the number of 

233 electrons, i.e., leads to p doping, which is accompanied by a lowering of the Fermi energy. 

234 The latter is tantamount to a shift of the graphene bands to lower binding energies with 

235 respect to the Fermi level, as found in the calculation and displayed in Fig. [5j 

236 By comparing Fig. [5] (c) and (d) the effect of doping for the case of graphene adsorbed 

237 top-fcc on Ni(lll) can by analyzed. It shows that the effect of doping is very similar to the 
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238 case of free-standing graphene. Again, the form of the graphene bands is somewhat changed 

239 and all graphene bands are shifted to lower binding energies. Now, however, this shift can 

240 not be simply explained by a lowering of the Fermi energy, because the latter is determined 

241 by the substrate bands, that is the nickel bands. In order to understand why doping with 

242 boron leads to a shift of the graphene bands to lower binding energies also in the case of 

243 graphene adsorbed on nickel, we disregard the change of the form of the bands for a moment. 

244 If we assumed that nickel and graphene bands, as well as their alignment, remains completely 

245 unchanged upon doping then the number of electrons on the graphene would not change 

246 because the Fermi level determined by the semi-infinite substrate would not change. This, 

247 however, would mean that the graphene layer would be negatively charged, by one electron 

248 per (2x2) unit cell, because the replacement of boron by carbon reduces the positive charges 

249 of the nuclei of the graphene sheet by this magnitude. The electrons charging the graphene 

250 would come from the nickel substrate. This means the nickel surface would be charged by the 

251 same magnitude as the graphene sheet, which would lead to a dipole layer. The electrostatic 

252 potential of such a dipole layer can be easily calculated as the potential of a plate capacitor 

253 and in the considered case of a charge of one electron per (2x2) unit cell would amount 

254 to a potential step of 17.5 eV between the nickel surface and the graphene layer (see, e.g., 

255 Ref. HI]). Raising of the graphene bands by 17.5 eV, of course, is completely unphysical and, 

256 indeed would result in an almost complete depopulation of the graphene valence bands and 

257 in a positive changing much higher than the initial negative charging of the graphene layer 

258 and, furthermore, would be in contradiction to the assumption that the alignment between 

259 nickel and graphene bands remains unchanged upon doping of graphene with boron. In order 

260 to avoid such a contradiction we now allow this level alignment to change in the course of 

261 doping. Then, the following can happen if we still disregard changes of the form of bands 

262 upon doping: Only a small amount of charge is transferred from the nickel surface to the 

263 doped graphene sheet. The potential step of this dipole layer then shifts the graphene bands 

264 to lower binding energies similarly as in the case of free-standing graphene. Indeed, shifts 

265 of the graphene bands of 1.75 eV, which is roughly the amount observed in free-standing 

266 graphene, would require only a transfer of 0.1 electron from nickel to graphene per (2x2) unit 

267 cell. Actually, in the considered scenario, the shift of the graphene bands to lower binding 

268 energies with respect to the Fermi level upon doping has to be somewhat smaller in the case 

269 of graphene adsorbed on nickel than in the case of free-standing graphene because otherwise 
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270 the graphene would be neutral again and no shift would occur. Therefore, one would expect 

271 a charging of the graphene layer upon doping by somewhat less than 0.1 electron per (2x2) 

272 unit cell. This, indeed, is found in our calculations. The charge of the pristine graphene on 

273 Ni (111), according to a Bader charge analysis, is 0.352 electrons per (2x2) unit cell, while 

274 that of the doped graphene on Ni(lll) is 0.424 electrons per (2x2) unit cell. This means 

275 doping leads to an increase of the negative charge of 0.072 electrons per (2x2) unit cell. 

276 In summary we can explain the shift of graphene bands to lower binding energies with 

277 respect to the Fermi level in the cases of free-standing and of graphene adsorbed top-fcc on 

278 Ni(lll) as follows: In the case of free-standing graphene the Fermi energy is lowered because 

279 of the smaller number of electrons due to doping, in the case of graphene adsorbed top-fcc 
2bo on Ni(lll) a small charge transfer from the nickel substrate to the graphene sheet occurs 
2bi which leads to dipole layers accompanied by a potential step between the nickel surface 

282 and the graphene sheet that is responsible for the shift of the graphene bands. This means 

283 doping has a similar effect for free-standing graphene and graphene adsorbed on nickel but 

284 the explanation for the effect is different in the two cases. 

285 Comparing Fig. [5] (a) and (c) and Fig. [5] (b) and (d) gives information about the effect 

286 of adsorbing pristine or doped graphene, respectively, top-fcc on nickel. In both cases a 

287 lowering of the graphene bands to stronger binding energies is found that is accompanied 

288 by n doping of the graphene sheet. According to Ref. HH (see also Refs. I3"8TI4T)|) . the reason 

289 for this band shift and the n doping is due to the surface dipole layer of the nickel substrate 

290 and a rearrangement of nickel surface charge known as pillow effectP^^. 

291 Considering the joint influence of adsorption and of boron doping together, Fig. [5] (a) vs. 

292 (d), adsorption leads to a shift of the graphene bands to higher binding energies, whereas 

293 doping has a reverse effect. For the considered case of a doping of 0.25 ML these super- 

294 imposed effects yield n and tt* bands that are located above the majority of the substrate 

295 bands around K. This results in a Dirac point that is estimated, as midpoint between tt and 

296 tt*, to be located 0.7 eV above the Fermi level and that is opened by 0.3 eV. 

297 D. Influence of the Boron Concentration 

298 In Figs. [2] (a)-(c) typical ARPES data of the preparation of boron-doped graphene on 

299 Ni(lll), with boron contents ranging from 0.045 to 0.28 ML, are shown. The band structure 
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300 of graphene, especially the tt band and its typical dispersion leading upward from the T point 

301 to the M and the K point is visible in all three cases. The main difference between the three 

302 cases is the energetic position of the minimum of the tt band at the F point with respect 

303 to the Fermi level. Fig. [3] shows that this position is measured to decrease from 10.1 eV for 

304 pristine graphene to 8.93 eV for maximally doped graphene (doping of 0.28 ML). 

305 The energy shift from the ARPES measurements is in good agreement with the shift 

306 observed in corresponding DFT band structure calculations (the underestimation of the ab- 

307 solute binding energies in the DFT band structures by about 0.2-0.4 eV is typically and 

308 reflects the inaccuracy of binding energies from DFT calculations) . This confirms the agree- 

309 ment between our explanations derived from calculations and our experimental results, i.e., 

310 it confirms that the boron-doped graphene was predominantly synthesized in the substi- 
3n tutional doping arrangement. The observed shift of the tt bands towards smaller binding 

312 energies is in line with the findings for nitrogen doped graphene, where the tt band is shifted 

313 towards higher binding energies for the substitutional doping geometry. In the band struc- 

314 ture calculations of free-standing graphene, we find a similar shift (not displayed) of 1.45 eV 

315 of the position of the tt band at the T point, with respect to the Fermi level, for the case of 

316 a boron coverage of 0.25 ML relative to pristine graphene. In the free-standing doped cases 

317 the energetic position of the Dirac point is estimated as midpoint between the graphene 

318 7i and tt* band at K, due to the absence of nickel bands. The bandwidth is estimated as 

319 difference between Ed and the energy of the tt band at T, which reveals a small increase of 

320 the width of the tt band of 0.26 eV. This shows, that the observed changes of the tt band 

321 at r are due to a shift of all graphene bands relative to the Fermi level and should not be 

322 mistaken with a decreased band dispersion. 



323 IV. CONCLUSION 

324 We have shown that boron-doped graphene can be grown on a Ni(lll) surface in a CVD 

325 process using the boron-containing precursor triethylborane or by segregation from a boron 

326 rich nickel crystal. Doping by boron leads to a strong shift of the graphene valence bands 

327 to lower binding energies. Additionally, the graphene tt band becomes more diffuse in the 

328 case of high boron-doping. Our results show a facile way of doping graphene with boron 

329 accompanied by a tuning of the band energies, up to shifts of 1.2 eV. It is found that the 

15 



330 effect of the adsorption of graphene on Ni(lll), which leads to a shift of the graphene bands 

331 to higher binding energies, is counteracted by introducing boron, resulting in an opened 

332 Dirac point in the unoccupied states for high boron coverages. Furthermore, our DFT 

333 calculations showed that boron-doped graphene prefers to adsorb in the top-fcc geometry, 

334 due to the strong preference of boron to adsorb in fcc-hollow sites. Due to the smaller 

335 bonding distance of boron compared to carbon, we predict a bending of the graphene layer 

336 in the case of low boron concentrations. 
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